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(54) Bezeichnung: VERFAHREN UND VORRICHTUNG ZUR BILDREKONSTRUKTION EINES RAUMVOLUMENS 



(57) Abstract: The invention relates to a method and a de- 
vice for reconstructing an image in a spatial volume on the 
basis of acquired projections, wherein during a reconstruc- 
tion step each acquired projection or each region of an ac- 
quired projection is fed once into a data processing system 
from a memory for the acquired projections and the inten- 
sity of a voxel of the reconstructed spatial volume image is 
updated during the reconstruction step for each voxel rele- 
vant projection or region. 

(57) Zusammenfassung: Die Erfindung betrifft ein 
Verfahren und einen Vorrichtung zur Bildrekonstruktion 
eines Raumvolumens aus akquirierten Projektionen, bei 
welchem wahrend eines Rekonstruktionsschrittes jede 
akquirierte Projektion bzw. jede Region einer akquirierten 
Projektion einmal einer Datenverarbeitung aus einem 
Speicher fur die akquirierten Projektionen zugefuhrt 
und die Intensitat eines Voxels des rekonstruierten 
Raumvolumenbildes wahrend des Rekonstruktionsschrittes 
je fur das Voxel relevanter Projektion bzw. Region 
aktualisiert wird. 
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VERFAHREN UND VORRI CHTUNG ZUR BILDREKONSTRUKTION EINES RAUMVO LUMENS 



Die Erfindung betrifft ein Verfahren und eine Vorrichtung zur Bildrekon- 
struktion eines Raumvolumens aus akquirierten Projektionen. 

5 Die Generierung von dreidimensionalen Daten aus Projektionen ist far 
den medizinisehen Bereich beispielsweise in Computertomographiesyste- 
men (CTs) verwirklicht. Fur die Bildkonstruktion aus den hierbei aufge- 
nommenen Projektionen werden in der Regel spezielle Multiprozessorsys- 
teme eingesetzt, um akzeptable Rekonstruktionszeiten von wenigen Se- 
10 kunden pro rekonstruierter Schicht zu erzielen, wobei in konventionellen 
CT-Geraten eine Aufhahme Schicht fur Schicht erfolgt, wahrend neuere 
CT-Gerate durch mehrzeilige Detektoren in der Lage sind, bis zu vier 
Schichten gleichzeitig zu scannen. 

Inzwischen besteht grofies Interesse, auch konventionelle Rontgensysteme 
15 fur die tomographische Bildgebung einzusetzen. Zum Einen kann auf- 
grund der Kegelstrahlgeometrie und der Verwendung ausgedehnter, 
zweidimensionaler Detektoren ein ganzes Volumen auf einmal gescannt 
werden. Zum anderen sind die Kosten fur konventionelle Rontgensysteme 
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geringer, insbesondere da nicht extra ein CT angeschafft werden muss, 
und die hierbei erreichbaren Auflosungen in der Regel holier. 

Fur die Bildrekonstruktion aus hierbei akquirierten Projektionen wird in 
der Regel ein Riickprojektionsalgorithmus genutzt. Hierbei wird fur jedes 
5 Projektionsbild die in den Projektionen enthaltene Grauwertinformationen 
gleichmaBig langst eines Strahl vom Projektionspixel zur Strahlenquelle 
verteilt und mit einem geometrischen Abstandsfaktor gewichtet. Unter 
anderem finden derartige Ruckprojektionen im Zusammenhang mit einer 
gefilterten Projektion Anwendung, bei welcher die Projektionsdaten zu- 

10 nachst ortsabhangig gewichtet und einer Filterung unter wor fen und dann 
in einen mit Null initialisierten Volumendatensatz riickprojiziert werden. 
Daruber hinaus kann die Riickproduktion auch iterative Verfahren umfas- 
sen, bei welchen in jedem Schritt eine Projektion- sowie eine Ruckprojek- 
tionsoperation enthalten sind. Nach Beendigung der Bildrekonstruktion 

15 konnen die Volumendaten aus einem Speicher fur das rekonstruierte 
Raumvolumenbild ausgelesen und je nach Anwendung sofort zweidimen- 
sional in Form von Schnitten oder auch dreidimensional visualisiert wer- 
den. 

Es hat sich herausgestellt, dass die hierbei benotigten Zeiten erheblich 
20 sind, beispielsweise kann ein Datenvolumen der GroBe 2563 mittels eines 
Feldkamp-Algorithmusses in ungefahr 15 Minuten rekonstruiert werden. 
Bei iterativen Rekonstruktionsverfahren verlangert sich diese Zeit unge- 
fahr urn die Anzahl der benotigten Iterationen, welche typischerweise in 
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der GroBenordnung einiger zehn Schritte liegt. Derartige Rechenzeiten 
sind jedoch von einer Echtzeit-Rekonstruktion weit entfernt und, insbe- 
sondere wenn die entsprechenden Einrichtungen im Zusammenspiel mit 
einer Positionskontrolle bzw. -uberwachung genutzt werden soli, nicht 
5 akzeptabel . 

Es bestehen Ansatze, die benotigten Rechenzeiten zu verkiirzen. 
Beispielsweise wird von T. Bortfeld in "optimized planning using physical 
objectives and constrains", Sem. In Rad. One, 9: 20-34 (1999) vorge- 
schlagen, bei einer Parallelstrahlgeometrie und einer gemeinsamen Dreh- 

10 achse, die Eigenschaft der Radontransformation auszunutzen, dass im Fu- 
rierraum die Projektion eine Ebene der Furiertransformierten des Volu- 
mens darstellt, so dass man die Rucktransformation durch Neuabtasten im 
Furierraum und Rucktransformation (Furier-Slice-Theorem) erhalt. Diese 
Moglichkeit ist jedoch auf lediglich die beschriebenen Geometrien be- 

15 schrankt. Dariiber hinaus besteht die Moglichkeit, die Rechengeschwin- 
digkeit zu erhohen, indem der Algorithmus parallelisiert wird. Nach Sa- 
saki T, Fukuda Y. in "reconstruction of 3-D X-ray computerized tomo- 
graphy images using a distributed memory multiprocessor system", trans- 
action of the information processing society of Japan, vol. 38, no. 9, 

20 Sept. 1997, Seiten 1681 bis 1693, und nach Dehner G, Herbert M. in 
"vector computing in CT image reconstruction-algorithm, data rate, com- 
pute power, parallel processing" in SPEEDUP, vol. 9, no, 2, Dec. 1995, 
Seiten 41 bis 47 (Schweiz) ist eine derartige Parallelisierung ohne Prob- 
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leme moglich. Es hat sich jedoch herausgestellt, dass bei einer derartigen 
Parallelisierung die Performance leider nahezu linear anwachst, da das 
Verfahren im Wesentlichen speicherbegrenzt ist. Die eigentliche Rechen- 
operation lauft in der Regel schneller als das Lesen der dafiir notwendigen 
5 Daten. Da die Rechengeschwindigkeit der Prozessoren immer groBer 
wird, die Speicherbandbreite (die Zahl der pro Zeit auslesbaren Daten) 
aber nur langsam zunimmt, wird dieses Verhaltnis ungiinstiger. Insofern 
ist durch die Parallelisierung theoretisch eine beliebige Performancestei- 
gerung moglich, was jedoch nur zu Lasten extremer Verteuerung der 

10 Hardware realisierbar ist. Eine Alternative besteht in der Nutzung 
spezieller Hardware fur die Rekonstruktion. Beispielsweise konnte von 
Ajakuijala J, Jaske UM, Sallinen S, Hehminen H, Laitinen J in 
"reconstruction of digital radiographs by texture mapping, ray casting and 
splatting", proceedings of the 18th Annual International Conference of the 

15 IEEE Engineering in Medicine and Biology Society" "Bridging Disci- 
plines for Biomedicine" Cat. No. 96CH36036) IEEE. Part vol. 2, 1997, 
Seiten 643 bis 645 vol. 2. New York, NY, USA dargelegt werden, dass 
mit Texture Mapping auf 3D-Graphikkarten (OpenGL-Implementierung) 
die CT-Rekonstruktion beschleunigt werden kann. Allerdings sind 

20 derartige Systeme nur bedingt zur Parallelverarbeitung geeignet und nicht 
zu skalieren. Weiterhin besteht nach Tresp V, Snell R, Gmitro AF in 
"videographic tomography II reconstruction with fan-beam projection 
data" IEEE transaction on medical imaging, vol. 13, no. 1, March 1994, 
Seiten 137 bis 143, USA die Moglichkeit, den rechenaufwendigen Teil 
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143 , USA die Moglichkeit, den rechenaufwendigen Teil optisch zu losen, 
was sich jedoch in der Praxis nicht durchsetzt. 

Es ist Aufgabe vorliegender Erfindung, ein Verfahren bzw. eine Vorrich- 
tung zur Bildrekonstruktion eines Raumvolumens aus akquirierten Projek- 
5 tionen bereit zu stellen, bei welchem unter minimalen Hardwareaufwand 
eine moglichst hohe Verarbeitungsleistung erreicht wird. 

Gelost wird diese Aufgabe einerseits durch ein Verfahren zur Bildrekon- 
struktion eines Raumvolumens aus akquirierten Projektionen, bei wel- 
chem wahrend eines Rekonstruktionsschrittes jede akquirierte Projektion 

10 einmal einer Datenverarbeitung aus einem Speicher fiir die akquirierten 
Projektionen zugefuhrt und die Intensitat eines Voxels des rekonstruierten 
Raumvolumenbildes wahrend des Rekonstruktionsschrittes je fiir das Vo- 
xel relevanter Projektion aktualisiert wird. Auf diese Weise ist es mog- 
lich, die Zahl der zu ladenden Projektionen zu minimieren und somit 

15 durch Reduktion der Anzahl der zeitaufwendigsten Speicheraktionen, 
namlich dem Einlesen von Projektionsdaten, die Gesamtzeit fur die Ver- 
fahrensdurchfiihrung auf ein Minimum zu reduzieren. 

Hierbei kommt es bei dem vorliegenden Verfahren jedoch nicht, wie bei- 
spielsweise bei dem Verfahren nach der DE 42 24 568 Al, darauf an, aus 
20 einem dreidimensionalen Datensatz ein zweidimensionales Schattierungs- 
bild zu entwickeln, sondern es soil aus den einzelnen Projektionen der 
dreidimensionale Datensatz ermittelt werden. 
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Je nach relativer Lage der Projektionsebenen zueinander kann es vorteil- 
haft sein, statt der jeweiligen Gesamtprojektionsebene lediglich Regionen 
einer Projektionsebene entsprechend zu behandeln. Dementsprechend wa- 
re es dann erforderlich, bei Veranderung der Lage der Voxel andere Re- 
5 gionen aus den jeweiligen Projektionen zu wahlen und fur den jeweiligen 
Rekonstruktionsschritt zu verwenden. Hierbei versteht es sich, dass je 
nach Lage dieser Regionen bestimmte Bereiche einer Projektionsebene 
mehrfach in dieser Weise behandelt werden, wobei erfindungsgemaB jede 
Region in seiner besonderen Form nur einmal geladen werden solL Insbe- 

10 sondere lasst sich hierdurch der gesamte Hardwarebedarf reduzieren, da 
fiir die Speicherung der in den Projektionen enthaltenen relevanten Infor- 
mation nur der hierfur unbedingt notwendige Speicherbedarf bereit ge- 
stellt werden muss. Insbesondere diese Vorgehensweise reduziert den 
Speicherbedarf fiir einen Zwischenspeicher bzw. einen Cache, in dem 

15 diese Informationen bereitgestellt werden konnen. 

Andererseits wird diese Aufgabe durch eine Vorrichtung zur Bildrekon- 
struktion eines Raumvolumens aus akquirierten Projektionen mit einem 
Speicher fiir die akquirierten Projektionen und einem Speicher fur das 
rekonstruierte Raumvolumenbild gelost, bei welcher diese Speicher durch 
20 eine Datenverarbeitung miteinander verkniipft sind und welche sich da- 
durch kennzeichnet, dass die Datenverarbeitung mindestens zwei Bearbei- 
tungspipelines umfasst, die einerseits je mit mindestens einem Speicherbe- 
reich eines Projektionscache fur eine Projektion oder eine Region einer 
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Projektion und andererseits mit mindestens einem Speicherbereich eines 
Voxelzwischenspeichers verbunden sind, wobei der Voxelzwischenspei- 
cher mit dem Speicher fur das rekonstruierte Raumvolumeribild und der 
Projektionscache mit dem Speicher fur die akquirierten Projektionen ver- 
5 knupft ist. 

Auf einer derartigen Vorrichtung ist das vorbeschriebene Verfahren bei- 
spielsweise durchfuhrbar , wobei die Verarbeitungsgeschwindigkeit schon 
deutlich tiber den Geschwindigkeiten herkommlicher Systeme liegt. Ande- 
rerseits ist unter Beibehaltung dieses Vorzuges auch eine andere Verfah- 
10 rensftihrung auf einer derartigen Anlage denkbar. 

Vorzugsweise werden somit mehrere Projektionen bzw. Regionen dieser 
Projektionen parallel verarbeitet, so dass diese lediglich einmal geladen 
werden miissen, wobei insbesondere bei einer derartigen Anordnung fur 
eine einmal geladenen Projektion bzw. Region einer Projektion alle rele- 
15 vanten Voxel entsprechend der Rekonstruktionsvorgaben bearbeitet wer- 
den. 

Insbesondere kann in einem Unterschritt der Rekonstruktion ein Subwiir- 
fel aus dem Voxelraum bearbeitet werden, wobei aus Effektivitats- und 
somit Geschwindigkeitsgriinden eine derartiger Subwurfel in dem Vo- 
20 xelzwischenspeicher abgelegt wird. Durch die Ablage des Subwiirfels in 
den Zwischenspeicher stehen diese Intensitatsdaten fur die Auswertung in 
dem Zwischenspeicher zur Verfugung und es ist nicht notig, die jeweili- 
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gen Daten aus dem eigentlichen Speicher fur das rekonstruierte Raumvo- 
lumenbild zu entnehmen, was wegen dessen GroBe normalerweise sehr 
zeitaufwendig ist. Fiir das Fiillen des Zwischenspeichers kann bei einem 
derartig gefiihrteii Verfahren ein gesonderter Algorithmus bzw. eine ge- 
5 sonderte Hardwarestruktur genutzt werden. 

Insbesondere im Zusammenspiel mit mehreren Bearbeitungspipelines 
bzw. der vorbeschriebenen Parallelrechnerstruktur - aber auch unabhan- 
gig hiervon - kann eine Bearbeitung des gesamten Raumvolumenbildes in 
Form geeignet gewahlter Subwiirfel auch unabhangig von der ubrigen 
10 Verfahrensfuhrung vorteilhaft die Verarbeitungsgeschwindigkeit erhohen. 
Dieses insbesondere dann, wenn der entsprechende Voxelzwischenspei- 
cher als Schieberegister ausgebildet ist. 

Wie bereits vorstehend angedeutet, kann es vorteilhaft sein, in einem Re- 
konstruktionsunterschritt zunachst eine Auswahl von Voxeln zu Bearbei- 

15 ten, deren Projektionen auf einer gleichen Region einer Projektionsebene 
abgebildet werden. Hierdurch lassen sich die in dieser Region vorhande- 
nen Projektionsinformationen in einem verhaltnismaBig kleinem Speicher, 
insbesondere einem Cache zur Verfugung stellen, wobei durch einen Ca- 
che, unabhangig von der ubrigen erfindungsgemaften Verfahrensfuhrung 

20 bzw. unabhangig von den ubrigen Merkmalen der beschriebenen Vorrich- 
tung, die Verarbeitungszeit verringert wird, da ein derartiger Cache we- 
sentlich geringere Zugriffzeiten fiir eine mit diesem verbundene Rechner- 
einheit aufweist. 
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Die entsprechende Region kann beispielsweise dadurch ermittelt werden, 
dass die jeweiligen Voxel jeweils auf die Projektionsebene projiziert und 
die jeweilige, liber deckte Flache genutzt wird. 

Vorzugsweise ist diese Auswahl von Voxeln ein vorbeschriebener Sub- 
5 wiirfel, wodurch sich die Vorteile der vorbeschriebenen Losungsansatze 
summieren. Es versteht sich, dass die Subwiirfel nicht zwingend eine 
Wurfelform aufweisen miissen, vielmehr kann jede Voxelmenge, vor- 
zugsweise jede einfach zusammenhangende Voxelmenge zur Anwendung 
kommen, wenn hierdurch die Zahl der Projektionen bzw. Regionen, die 
10 hintereinander fiir eine Folge bestimmter Rekonstruktionsschritte benotigt 
wird, reduziert werden kann. 

Die Verwendung eines Voxelzwischenspeichers ist auch unabhangig von 
einem Pr oj ektionscache bzw. einer Parallelrechnerstruktur vorteilhaft, da 
sich durch einen derartigen Zwischenspeicher die Zugriffszeiten der Re- 
15 cheneinheit verkurzen las sen, da ein groBer Hauptspeicher, wie er fur die 
Auftiahme der gesamten, rekonstruierten Bildinformation benotigt wird, 
wesentlich langsamere Zugriffsgeschwindigkeiten aufweist. 

Um die Bearbeitungszeit weiter zu reduzieren, kann ein zweiter Vo- 
xelzwischenspeicher vorgesehen sein, der alternativ zu dem ersten Vo- 
20 xelzwischenspeicher mit der Bearbeitungspipeline verbunden werden und 
unabhangig von dem ersten Voxelzwischenspeicher mit dem Speicher fur 
das rekonstruierte Raumvolumenbild Daten austauschen kann. Insofern 
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kann einer der beiden Voxelzwischenspeicher mit dem Speicher far das 
rekonstruierte Raumvolumenbild Daten austauschen, wahrend der andere 
Voxelzwischenspeicher fur die Rechenoperationen genutzt wird. 1st der 
Rechenvorgang beendet, so konnen durch einen einfachen Umschalter 
5 diese Verbindungen getauscht werden. Auf diese Weise werden Totzeiten 
durch den Datenaustausch zwischen Voxelzwischenspeicher und Speicher 
fur das rekonstruierte Raumvolumenbild vermieden. 

Des weiteren schlagt die Erfindung eine Vorrichtung zur Bildrekonstruk- 
tion eines Raumvolumens aus akquirierten Projektionen mit einem Spei- 

10 cher fiir die akquirierten Projektionen und einem Speicher fur das rekon- 
struierte Raumvolumenbild, die iiber eine Datenverarbeitung miteinander 
verkmipft sind, vor, bei welcher die Speicherbandbreite unter der Verar- 
beitungsleistung des Gesamtsystems liegt. Eine derartig ausgerichtete 
Vorrichtung ist in der Lage, schneller zu arbeiten, als dieses die Speicher 

15 ermoglichen und bedingt somit, dass die Hardware optimal genutzt wird. 
Vorzugswiese werden Speicherbandbreite und Verarbeitungsleistung in 
Voxel/Sekunde verglichen, wobei hier auch andere Kriterien, die einen 
Vergleich zwischen der Leistung der Speicher fur das rekonstruierte 
Raumvolumenbild und fiir die akquirierten Projektionen mit der Verarbei- 

20 tungsleistungen erlauben, moglich sind. 

Es versteht sich, dass in vorliegendem Zusammenhand der Begriff 
"Raumvolumenbild" jede Darstellung umfasst, in der in den Projektionen 
enthaltene Information dreidimensional riickermittelt und gespeichert ist. 



WO 02/061686 



11 



PCT/DE02/00266 



Insbesondere kann es sich hierbei urn eine Intensitatsverteilung in einem 
Voxelraum handeln. Selbiges gilt auch fur die "Projektionen". 

Weitere Eigenschaften, Ziele imd Vorteile der Erfindung werden anhand 
nachfolgender Beschreibung anliegender Zeichnung erlautert, in welcher 
beispielhaft erfindungsgemaBe Verfahrensschritte sowie eine Rechnerar- 
chitektur erlautert ist. In der Zeichnung zeigen: 

Fig. 1 eine schematische Darstellung einer Rontgenanlage; 

Fig. 2 die Rontgenanlage nach Fig. 2 im Schnitt; 

Fig. 3 eine schematische Rechnerarchitektur fur die Ruckpr oj ektion ; 

Fig. 4 eine weitere Rechnerarchitektur fur die Ruckproj ektion; 

Fig. 5 eine schematische Darstellung der Rechnerstruktur nach Fi- 

gur 3 mit der Verkniipfung zwischen Volumen und Projekti- 
on; 

Fig. 6 einen Verfahrensablauf mit der Rechnerarchitektur nach den 

Figuren 3 bis 5; 

Fig. 7 die Auswahl eines Subwiirfels; 

Fig. 8 die Auswahl geeigneter Regionen fiir eine Anzahl mehrere 

Subwurfel bzw. Voxel; 

Fig. 9 mehrere Projektionsebenen fur einen Subwurfel; 
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Fig. 10 eine mogliche Anordnung von Projektionen und zu rekon- 

struierendem Bildraum; 

Fig. 11 die Auswahl einer Voxelscheibe, einer in der Voxelscheibe 

gewahlten Voxelzeile bzw. eines in der Voxelzeile gwahlten 
5 Voxel wiirfels ; 

Fig. 12 eine Rechnerstruktur nach dem Stand der Technik mit mehre- 

ren parallelen Projektionsspeichern; und 

Fig. 13 eine Rechnerstruktur nach dem Stand der Technik mit mehre- 

ren parallelen Speichern fur die Intensitatsinformationen. 

10 Bei der in Fign. 1 und 2 schematisch dargestellten Rontgenanlage wird 
eine Person 1 mittels einer Strahlenquelle 2 durchstrahlt. Hierdurch kon- 
nen mit einem entsprechenden Detektor Projektionen 3 aufgenommen 
werden, die letztlich die Wechselwirkung des entsprechenden Strahlenke- 
gels 4 mit dem Korper der durchstrahlten Person widerspiegeln. Die 

15 Strahlenquelle 2 und der entsprechende Detektor sind um die Person ro- 
tierbar angeordnet, so dass verschiedene Projektionsrichtungen aufge- 
nommen werden konnen. Es versteht sich, dass start einer derartigen Ein- 
richtung auch andere Einrichtungen, bei denen eine Riickprojektion not- 
wendig ist, zur Anwendung kommen konnen. Insbesondere konnen auch 

20 Gegenstande entsprechend untersucht werden. 
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Die ermittelten Projektionen 5 werden in einem entsprechenden Speicher 
6 fur die akquirierten Projektionen abgelegt Aus diesen soli ein Raumvo- 
lumenbild ermittelt werden, welches in einem Speicher 7 fur das rekon- 
struierte Raumvolumenbild abgelegt wird (siehe Figuren 1 bis 5). 

5 Die beiden Speicher 6 und 7 sind uber eine Datenverarbeitung 8 mitein- 
ander verkniipft, wobei bei den in Figuren 3 und 4 dargestellten Ausfuh- 
rungsbeispielen die Datenverarbeitung 8 einen Projektionscache 9 um- 
fasst. In diesen Projektionscache 9 konnen nach Bedarf Daten aus dem 
Projektionsspeicher 6 abgelegt werden. Der Projektionscache 9 weist ein- 

10 zelne Speichersegmente 10 auf, hi welche jeweils ubermittelte Pixeldaten 
einer Projektionsebene abgelegt werden. Je Speicher einheit 10 ist eine 
Hardware- beziehungsweise Bearbeitungspipeline 11 vorgesehen, welche 
einer Zelle 12 eines als Schieberegister ausgebildeten Voxelzwischenspei- 
cher 13 zugeordnet ist. Die Hardwarepipeline 11 liest fur das in dem 

15 Speicher 12 abgelegte Volumenelement die notwendigen Projektionspixel 
aus dem Projektionsspeicher 9. Darin berechnet sie den Beitrag fur das 
Volumenelement, der zu dem bisherigen Beitrag des Volumenelementes 
zuaddiert wird. Hierbei werden im vorliegenden Ausfuhrungsbeispiel zur 
Ruckprojektion die Intensitaten langs der Strahlen „ Quelle- 

20 Projektionspixel" gleichmaBig iiber die betreffenden Voxel verteilt, wobei 
- je nach konkreter Ausfiihrungsform - noch eine Funktion vorgesehen 
sein kann, welche eine geometrische Abschwachung berucksichtig. 
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Hierbei muss fiir jedes Voxel ausgehend von jeder Projektion berechnet 
werden, auf welchen Punkt das Zentrum des Voxels abgebildet wird. Die 
Intensitat des Punktes in der Projektionsebene 5 wird in der Regel durch 
bilineare Interpolation der benachbarten Pixelintensitaten berechnet. Der 
5 ermittelte Wert wird dann mit dem inversen Quadrat des Abstandes 
„ Quelle- Voxel" multipliziert und zum bisherigen Beitrag im Voxel ad- 
diert. 

Da alle Hardwarepipelines 11 synchron arbeiten, werden sie zur gleichen 
Zeit fertig. Daraufhin werden die Ergebnisse in dem Schieberegister 13 
10 ein Speicherelement 12 weiter geschoben, sodass fur jedes Volumenele- 
ment der Beitrag der nachsten Projektionsebene bestimmt wird. 

Sind die Beitrage aller Projektionsebenen berechnet und akummuliert 
worden, so werden die jeweiligen Daten wieder in den Speicher 7 zu- 
riickges chrieben . 

15 Um die Zeit fiir das Rixckspeichern sowie das neue Abspeichern der je- 
weiligen Daten von dem Voxelzwischenspeicher 13 zu dem Speicher 7 
nicht ungenutzt verstreichen zu lassen, ist bei der Ausfuhrungsform nach 
Figur 4 ein zweiter Voxelzwischenspeicher 14 vorgesehen, der wahrend 
der Umspeichervorgange des ersten Voxelzwischenspeicher s 13 mit den 

20 Hardwarepipelines 11 verbunden wird. Auf diese Weise werden die Pa- 
rallelprozessoren der Hardwarepipelines 1 1 optimal genutzt. Wahrend der 
Berechnung der in dem Zwischenspeicher 14 gespeicherten Voxel wird 
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der Voxelzwischenspeicher 13 entsprechend geleert und neu geladen, so- 
das s er dann fur weitere Berechnungen zur Verfugung steht, sowie die 
Berechnungen fur die Daten in dem Voxelspeicher 14 abgeschlossen wer- 
den. Daraufhin wird der Voxelspeicher 13 wieder mit den Hardwarepipe- 
5 lines 11 verbunden, wahrend fur den Voxelspeicher 14 der Datenaus- 
tausch mit dem Speicher 7 vorgenommen wird. 

Zwar ist es denkbar, dass jeweils samtliche Daten einer Projektionsebene 
5 in den Projektionscache 9 eingelesen werden. Da jedoch das Volumen 
des Voxelzwischenspeichers 13 begrenzt ist, wurde eine groBe Zahl re- 

10 dundanter Daten geladen werden, weil je Projektion far ein bestimmtes 
Voxel nur sehr kleine Regionen relevant werden konnen. Insofern ist es 
vorteilhaft, nur derartige Regionen, wie sie beispielsweise in Figur 5 mit 
der Ziffer 15 bezeichnet sind, zu laden. In vorliegendem Zusammenhang 
wird unter dem Begriff „ Region" eine verhaltnismaBig geringe Menge an, 

15 vorzugsweise einfach zusammenhangenden, Pixeln verstanden, wobei die 
GroBe der Menge derart gewahlt ist, dass diese Pixelmenge ohne weiteres 
in den Projektionscache 9 geladen werden kann. 

Prinzipiell konnen in den Projektionscache 9 und in dem Voxelzwi- 
schenspeicher 13 bzw. 14 beliebige Voxel bzw. Projektionen 5 oder Re- 
20 gionen 15 geladen werden. Vorzugsweise werden jedoch die in den Pro- 
jektionscache 9 eingelesenen Regionen 15 miteinander korreliert. Die 
Kor relation kann derart gewahlt werden, dass aus dem zu rekonstruier en- 
den Bildraum 16 ein Subwiirfel 17 gewahlt wird, dessen Zahl an Voxeln 
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vorzugsweise der Zahl der Speicherelemente 12 in dem Voxelzwi- 
schenspeicher 13 bzw. 14 entspricht. Bespielhaft ist in Figur 7 ein Sub- 
wiirfel 17 mit vier Voxeln Kantenlange dargestellt, der in einem Vo- 
xelzwischenspeicher 13 mit 64 Speicherelementen abgelegt werden kann. 

5 In den Projektionscache 9 werden dann alle die Regionen 15 geladen, 
welche die jeweiligen relevanten Bilddaten korrespondierend zu der je- 
weiligen Projektionsrichtung enthalten, wie dieses beispielhaft in Figuren 
5 und 9 dargestellt ist. Sowie die entsprechenden Projektionsdaten geladen 
sind, konnen die jeweiligen Berechnungen durchgefiihrt werden, wobei 

10 durch das Schieberegister parallel samtliche Voxel des Subwiirfels 17 fur 
die jeweilige in dem Projektionscache 10 geladenen Regionen 15 berech- 
net werden, Zwar sind idealer Weise in dem Projektionscache 9 samtliche 
Projektionsrichtungen abgelegt. Es kann jedoch eine Kosten-Nutzen- 
Abwagung vorgenommen werden, nachdem lediglich ein Teil der not- 

15 wendigen Projektionen 5 in dem Projektionscache 9 abgelegt ist und zwi- 
schenzeitlich ein entsprechender Datenaustausch vorgesehen ist. 

Der somit durchgefiihrte Verfahrensablauf ist in Figur 6 dargestellt. Zu- 
nachst wird ermittelt, welche Regionen 15 der Projektionsebene 5 zu ei- 
nem entsprechenden Subwurfel 17 beitragen. Dieses kann bei der Pro- 
20 gramminstallation bzw. bei der Konzipierung der Hardwarestruktur be- 
reits erfolgen. AnschlieBend werden jeweils far einen Subwurfel in einer 
inner en Schleife samtliche Beitrage der Projektionen fur ein entsprechen- 
des Voxel ermittelt, wobei der Subwurfel 17 mittels des Schieberegisters 
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13 abgearbeitet wird. AnschlieBend wird ein weiterer Subwiirfel 17 ge- 
wahlt, wobei nach Moglichkeit die hierbei beriicksichtigten Regionen bei- 
behalten werden und nur selektiv einzelne Regionen 15 ausgetauscht wer- 
den brauchen. Gegebenfalls wird aller dings der Satz an Regionen 15 
5 komplett ausgetauscht. 

Da bei dem vorliegenden Ausfiihrungsbeispiel davon ausgegangen wird, 
dass die Zahl der Projektionen 5 die Zahl der Speicher 10 des Projekti- 
onscache 9 ubersteigt, wird, wenn samtliche Subwiirfel 17 des Volumens 
16 abgearbeitet sind, ein kompletter neuer Satz Projektionsrichtungen 
10 bzw. Projektionen 5 gewahlt, und es werden die Beitrage dieser Projekti- 
onen entsprechend bestimmt. 

In einer konkreten Ausfuhrungsvariante basiert das Rekonstruktions- 
verfahren im wesentlichen auf zwei Schritten: zuerst werden die Daten 
gefiltert, dann die Ruckprojektion durchgefuhrt. Ohne Beschrankung der 

15 Allgemeinheit kann man dabei annehmen, dass sich die Rontgenanlage um 
eine Achse gedreht hat, die parallel zu der Normalen einer der Seitenfla- 
chen des Volumenwiirfels bzw. des zu rekonstruierenden Bildraums 16 
liegt (vergl. Fig. 10). Im Ruckprojektionsschritt konnen dann Voxelschei- 
ben 17A der Dicke N parallel zu dieser Seitenflache getrennt betrachtet 

20 und aus entsprechenden Projektionszeilen 5A (exemplarisch in Fig. 10 
beziffert) der Projektionen 5 rekonstruiert werden. Aus diesen Scheiben 
konnen entsprechend Fig. 11 Voxelwurfel 17 der GroBe NxNxN gelesen 
werden, wobei dieses vorzugsweise iterativ nach Voxelzeilen 17B erfolgt. 
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Der Beitrag eines Stapels von NxNxN Projektionen (Projektionsblock) auf 
jedes der Voxel dieses ausgelesenen Wurfels wird berechnet und auf die 
aktuellen Werte der Voxel addiert. 

Die Projektionen 5 konnen dann jeweils in ihren fur eine Voxelscheibe 
5 17 A relevanten Zeilen 5 A zur Berechnung der einzelnen Voxel genutzt 
werden. Insbesondere kann, wenn die Teile der Projektionen 5, welche 
die jeweilige Voxelwiirfel-Scheibe 17A betreffen, fur alle Projektionen 5 
geladen sind, der Beitrag zu der gesamten Voxelwiirfel-Scheibe 17A be- 
rechnet werden. Daher enthalt vorzugsweise jede Pipeline 11 genugend 

10 Speicher fur diese Daten. Aufgrund dieser geeigneten Anordnung von 
Projektion zu Volumen kann man davon ausgehen, dass es sich um meh- 
rere benachbarte Zeilen in dem Projektionsbild handelt. Damit das Ver- 
fahren noch effizienter arbeitet, werden vorzugsweise die einzelnen Zei- 
len wahrend der Berechnung nachgeladen (2-way -memory) , um Zugriffs- 

15 konflikte zu vermeiden, bzw. mehrere Speicherbanke genutzt. 

Die GroBe des Speichers betragt bei dieser Ausfuhrungsvariante vor- 
zugsweise ca. Voxel wtirfelbreite * Seitenlange des Projektionsbildes (ty- 
pischerweise 512 oder 1024) Werte (typischerweise 16 Bit). Die Voxel- 
wtirfelbreite ist typischerweise 4 ( da 64 Pipelines 4x4x4 entsprechen) , 
20 d.h. es ist ffir jede Pipeline ein Puffer der GroBe 4x2 Byte x 1024 = 8 
kByte notwendig, insgesamt also 512 kByte Pufferspeicher. Fiir jede der 
64 Pipelinestufen werden vorzugsweise 8 Multiplizierer genutzt, d.h. 512 
Multiplizierer insgesamt. 



WO 02/061686 



19 



PCT/DE02/00266 



Der Berechnungsablauf sieht dann vorzugsweise wie folgt aus: 

Fiir jede Scheibe der Dicke N // Ruckprojektionsschritt 

Fur jeden Projektionsblock 

schreibe die Zeilen in das Projektionscache, die fiir die Schei- 
be benotigt werden 

Fur jeden Wiirfel der GroBe NxNxN der Scheibe 

lese den Wiirfel in das Schieberegister mit NxNxN Re- 
gister 

fuhre NxNxN mal aus 

schiebe die Voxel im Schieberegister um eine Stel- 
le weiter 

Parallel: fiir jedes Schieberegister 

berechne den Beitrag der dazugehorigen 
Projektionsebene zum Voxel im Register 

addiere diesen Beitrag auf das Voxel 

schreibe Wiirfel auf Volumenspeicher zuriick 

Insofern konnen insbesondere gleichzeitig Wiirfel und Zeilen von mehre- 
ren Projektionen gleichzeitig geladen gehalten und iiber den Schieberegis- 
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termechanismus jedes Voxel des Wurfels mit jeder der geladenen Projek- 
tionen zusammengebracht werden. Dieses hat den Vorteil, dass die be- 
rechnete Scheibe nach Durchlauf durch alle Projektionen fertig ist und 
fast alle verwendeten Zeilen der Projektionen nicht mehr geladen werden 
5 mussen. 

Beispielhaft sind in den Figuren 12 und 13 Ausfuhrungsbeispiele nach 
dem Stand der Technik dargestellt, wobei diese jeweils parallele Struktu- 
ren 108 bzw. 208 aufweisen, jedoch entweder die Projektionsebenen 105 
in mehr er en Speichern 106 oder aber die Voxelraume 207 in Parallelspei- 
10 chern abgelegt sind, wodurch die Kosten erheblich steigen. 

Vorliegende Erfindung kann beispielsweise in einer C-Arm- 
Angiographieanlage oder an eine Linearbeschleunigung in Verbindung 
mit electronic portal imagine device (EPID) umgesetzt werden. Die Re- 
konstruktion kann beispielsweise durch eine gefilterte Ruckprojektion 

15 bzw. durch imperative Verfahren erfolgen. Fur die gefilterte Ruckprojek- 
tion werden die Projektionsdaten zunachst ortsabhangig gewichtet und 
einer Filterung unterworfen. Bei einer Implementierung der Filterung im 
Frequenzraum unter Anwendung einer gut optimierten Software zur Fou- 
riertransformation kann der Schritt als relativ zeitunkritisch angesehen 

20 werden. Zur Ruckprojektion der gefilterten Profile werden diese sowie 
ein mit null initialisierter Volumendatensatz zunachst in den Speicher der 
entsprechenden Anlagenkarte geladen. Die Ruckprojektion erfolgt wie 
vorstehend beschrieben. Bei dem interativen Verfahren zur Bildrekon- 
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struktion ist in jedem Schritt eine Projektions- sowie eine Ruckprojekti- 
onsoperation enthalten. Zur Realisierung einer effizienten Projektionsope- 
ration (raytracing) sind bereits spezielle Verfahren und Architekturen be- 
kannt. Fur eine voxelbasierte Rtickposition kann das vorstehende be- 
5 schriebene Verfahren sowie die vorstehend beschriebene Vorrichtung 
zum Einsatz kommen. 

Durch die Verwendung der vorgeschriebenen Architektur sowie des vor- 
geschriebenen Verfahrens konnen grofte Datenmengen rechtzeitig verar- 
beitet und rekonstruiert werden, sodass mit dem beschriebenen System 
10 mittels Hochenergie-CT beispielsweise online die Positionierung eines 
Patienten in einem Linearbeschleuniger verifiziert werden kann. 
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Patentanspriiche : 

1. Verfahren zur Bildrekonstruktion eines Raumvolumens aus akqui- 
rierten Projektionen, bei welchem wahrend eines Rekonstruktions- 
schrittes jede akquirierte Prqjektion bzw. jede Region einer akqui- 
rierten Projektion einmal einer Datenverarbeitung (8) aus einem 
Speicher (6) fur die akquirierten Projektionen zugefiihrt und die In- 
tensity eines Voxels des rekonstruierten Raumvolumenbildes wah- 
rend des Rekonstruktionsschrittes je fur das Voxel relevanter Pro- 
jektion bzw. Region aktualisiert wird. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass mehrere 
Projektionen bzw. Regionen von Projektionen parallel verarbeitet 
werden. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass 
die Projektionen bzw. Regionen von Projektionen zur Verarbeitung 
in einen Projektionscache (9) geladen werden. 

4. Verfahren nach einem der Anspriiche 1 bis 3 , dadurch gekenn- 
zeichnet, dass in einem Unterschritt der Bildrekonstruktion ein 
Subwurfel (17) bearbeitet wird, dessen Intensitatsinformation in ei- 
nem Voxelzwischenspeicher (13) abgelegt und wahrend der Bear- 
beitung anhand von Daten der Projektionen aktualisiert wird. 
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5. Verfahren nach einem der Anspriiche 1 bis 3, dadurch gekenn- 
zeichnet, dass in einem Unterschritt der Bildrekonstruktion eine 
Auswahl (17) von Voxeln bearbeitet wird, der en Projektion auf ei- 
ner gleichen Region (15) einer Projektionsebene abgebildet werden. 

6. Vorrichtung zur Bildrekonstruktion eines Raumvolnmens (7) aus 
akquirierten Projektionen (5) mit einem Speicher (6) filr die akqui- 
rierten Projektionen und einem Speicher (7) fur das rekonstruierte 
Raumvolumenbild, die iiber eine Datenverarbeitung (8) miteinander 
verkmipft sind, dadurch gekennzeichnet, dass die Datenverarbei- 
tung (8) mindestens zwei Bearbeitungspipelines (11) umfasst, die 
einer seits je mit mindestens einem Speicherbereich (10) eines Pro- 
jektionscache (9) und anderer seits je mit mindestens einem Spei- 
cherbereich (12) eines Voxelzwischenspeichers (13) verbunden ist, 
wobei der Voxelzwischenspeicher (13) mit dem Speicher (7) fiir das 
rekonstruierte Raumvolumenbild und der Projektionscache (9) mit 
dem Speicher (6) fiir die akquirierten Projektionen verkniipft ist. 

7. Vorrichtung nach Anspruch 6, dadurch gekennzeichnet, dass der 
Voxelzwischenspeicher (13) als Schieberegister ausgestaltet ist. 

8. Vorrichtung nach Anspruch 6oder 7, gekennzeichnet durch wenigs- 
tens einen zweiten Voxelzwischenspeicher (14) der alternativ zu 
dem ersten Voxelzwischenspeicher (13) mit den Bearbeitungspipe- 
lines (11) verbunden werden und unabhangig von dem ersten Vo- 
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xelzwischenspeicher (13) mit dem Speicher (7) fur das rekonstruier- 
te Raumvolumenbild Daten austauschen kann. 

9. Vorrichtung zur Bildrekonstruktion eines Raumvolumens (7) aus 
akquirierten Projektionen (5) mit einem Speicher (6) fixr die akqui- 
5 rierten Projektionen und einem Speicher (7) fur das rekonstruierte 

Raumvolumenbild, die tiber eine Datenverarbeitung (8) miteinander 
verkmipft sind, dadurch gekennzeichnet, dass die Speicherbandbrei- 
te unter der Verarbeitungsleistung liegt. 
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